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I. INTRODUCTION 



Exact or approximate solutions of the relativistic/non-relativistic wave equations have re- 
ceived great attentions. So far the solutions are in general obtained for the case of constant 
mass or at most time-dependent mass [1,2]. The effective mass solutions have received much 
attentions recently. A quite general hermitian effective Hamiltonian is used to describe the 
non-relativistic systems, such description is applied to study the semiconductor nanostruc- 
tures [3] . Another interesting problem is that the correct form of the kinetic energy operator 
for such a Hamiltonian, since the momentum, and the mass operators are no longer commute 
in the case of position-dependent mass, which is related to the problem of ordering ambigu- 
ity [4]. There are some important problems related to the ordering ambiguity concept, such 
as the dependence of nuclear forces on the relative velocity of the two nucleons [5, 6], the 
impurities of crystals [7]. In addition, many authors have studied to propose some effective 
Hamiltonians for non-relativistic case taking into account the dependence of the mass on 
position [8]. 

There are many efforts about solving the Schrodinger equation for the case of position- 
dependent mass by using different methods or schemes for different potentials, such as expo- 
nential type potential [4], Natanzon potentials by using a group-theoretical method [9], so- 
lutions in the case of mappings of the Morse+oscillator+ Coulomb potential [10], hyperbolic- 
type potentials [11], Morse, and Coulomb potential with the position-dependent mass [12, 
13], PT-symmetric anharmonic oscillators [14], the Morse-like potential in the scheme of 
super symmetric quantum mechanics [15], Kratzer and Scarf II potentials [16], deformed 
Rosen-Morse, and Scarf potentials [17]. Many authors have been also solved the Klein- 
Gordon, and Dirac equation by taking a suitable mass distributions in one and/or three 
dimensional cases for different potentials, such as Coulomb potential [18], Lorentz scalar 
interactions [19], hyperbolic-type potentials [20], Morse potential [21], and Poschl-teller po- 
tential [22]. 

Here we intend to solve the Klein-Gordon equation within the framework of an approx- 
imation to the centrifugal potential term. We study the effect of the mass varying with 
position on the energy spectra, and the eigenfunctions of the vector, and scalar Hulthen 
potential [23], which is widely used in nuclear, particle physics, atomic physics, condensed 
matter, and chemical physics [24-26]. For our task, we use a general parametric form of 
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the Nikiforov-Uvarov (NU) method, which is based on turning of a second order differential 
equations to a hypergeometric type equation [27]. 

The organization of this work is as follows. In Section II, we give briefly the parametric 
generalization of the NU-method. In Section III, we give the energy eigenvalue equation, 
and corresponding eigenfunctions for the vector, scalar Hulthen potential for any ^-values 
in the position-dependent mass background. We obtain also the results for the case of the 
constant mass, and we summarize our concluding in Section IV. 



II. NIKIFOROV-UVAROV METHOD 



The Schrodinger equation can be transformed into a second order differential equation 
with the following form 



2 , ,d 2 ^(s) . w ,d^(s) . n 

where cr(s), a(s) are polynomials, at most, second degree, and f(s) is a first degree polyno- 
mial. In order to find a particular solution, we take the following form 



(2) 



We get from Eq. (1) 



.M^ + r( S )^§W) = 0, (3) 
where <f>(s) can be written in terms of Rodriguez formula 



R d n 

*•(») = ^5? [""W^"' (4) 



and the weight function p(s) satisfies 



da(s) a(s) dp(s) 
ds p(s) ds 



(5) 



The other factor of the solution is defined as 



1 1p(s) 7r(s) 



ip(s) ds a(s) 

In the method, the polynomial n(s), and the parameter k are defined as [27] 



(6) 



n(s) = \ [a'(s) - f(s)] ± {\[a'(s) - f {s)f - a(s) + ka(s)} V \ (7) 



and 



A = k + tt'(s) . (8) 

where A is a constant, and given in Eq. (3). Since square root in the polynomial 7r(s) in Eq. 
(7) must be a square then this defines the constant k. Replacing k into Eq. (7), we define 



T(s)=f(s) + 2n(s). (9) 

Since p(s) > and a(s) > 0, hence the derivative of r(s) should be negative [27], which 
leads to the choice of the solution. If A in Eq. (8) is 

A = K = -„t> - „ = ,1,2,... (10) 

the hypergeometric type equation has a particular solution with degree n. 

In order to explain the general parametric form of the NU method, let us take the general 
form of a Schrodinger-like equation including any potential 



[s(l - «3s)] 2 ^^ + [s(l - a 38 )( ai - «2*)]^^ + Kis 2 + &a - = 0. (11) 

When Eq. (11) is compared with Eq. (1), we get 



f(s) = «i - a 2 s ; <t(s) = s(l - a 3 s) ; a(s) = -^s 2 + &s - 6- ( 12 ) 
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Substituting these into Eq. (7), we get 



7r(s) = «4 + 0:5s ± y («6 — ka 3 )s 2 + (07 + fc)s + as , (13) 
where the parameters in the above equation are as follows 

«4 = |(1 - OLi) , oi b — \ (a 2 - 2o 3 ) , 

a6 = «5 + 6, , «7 = 2a 4 a 5 - £ 2 , (I 4 ) 
«8 = «1 + £3 • 

In NU-method, the function under square root must be the square of a polynomial, so 

h,2 = -("7 + 2o 3 o 8 ) ± 2 v /a 8 a 9 , (15) 

where 

ctg = 030:7 + OgOg + «6- (16) 

The function ir(s) becomes 

7r(s) = a 4 + a 5 s - [(y r ag + a 3y /a^)s - y/o^} . (17) 

for the k- value k = — (07 + 20:30:8) — 2^/osOg , where we have to say that the different fc's 
lead to the different 7r(s)'s. We also have from Eq. (9) 

t(s) = oti + 2o 4 - (o 2 - 2o 5 )s - 2 [(v/ag" + a 3y /a^)s - y/a^] . (18) 

Thus, we impose the following for satisfying the condition that the derivative of r(s) must 
be negative 



r '( s ) = -(02 - 2o 5 ) - 2( v /«^ + o 3v /a^) 

= -2o 3 - 2( v /q^ + o 3v /a^) < 0. 
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(19) 



From Eqs. (8), (9), (18), and (19), and equating Eq. (8) with the condition that A should 
satisfy given by Eq. (10), we obtain 

a 2 n - (2n + l)a 5 + (2n + l)(^/al^ + a 3 ^/a^) + n(n - l)a 3 

+ a 7 + 2a 3 a 8 + 2y/a 8 a 9 = 0. (20) 

which is the energy eigenvalue equation of a given potential. 

Now, let us look the eigenfunctions of the problem with any potential. We obtain the 
second part of the solution from Eq. (4) 

Ms) = P^ 1 ^ aW ~ 1 \l-2a 3 s), (21) 
by using the explicit form of the weight function obtained from Eq. (5) 

p(s) = s^-^l - ass)^-" 10 ' 1 , (22) 

where 

ocw = o-i + 2a 4 + 2y/ag~ ; an = a 2 — 2a 5 + 2(^fau ) + a^y/as) . (23) 
and P( a ^\l — 2a 3 s) are Jacobi polynomials. From Eq. (6), one gets 

7p(s) = s ai2 (l-a 3 sy ai2 -^t , (24) 
then the general solution ^f(s) = ip(s)(f)(s) becomes 

= s «»(i _ assr^-^P^ 10 ' 1 '^'^^ - 2a 3 s) , (25) 

where 

«i2 = «4 + Vas ; «i3 = «5 - (v^9 + )• (26) 
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III. BOUND-STATE SOLUTIONS 



The Klein-Gordon equation for a particle with mass m with vector V v (r), and scalar V s (r) 
potentials is (h — c — 1) 

{ — V 2 — [E 2 — m 2 (r)] + 2[m(r)V s (r) + EV v {r)\ + [V 2 (r) - V?(r)]}*(r, 0, 0) = , (27) 
Using \I/(r, 9, 4>) = r^ 1 (f)(r)Y em (9, 0), we have the radial part o the equation 

f^l + {[E 2 -m 2 (r)}-2[m(r)V s (r) + EV v (r)} 

+ [V s 2 (r) - V 2 {r)] - }0(r) = . (28) 

where Yi m (9,<f)) is spherical harmonics, and £ is the angular momentum quantum number. 
In order to solve the Eq. (30), we prefer to use the following mass function 



m(r) =m + - _ ^_ r/rQ , (29) 

where m , mi are two arbitrary, positive constants. We have to use an approximation, given 
by 1/r 2 rs e r / ro / (e r / ro — l) 2 ^, to the centrifugal term, since the radial equation has no 
analytical solutions for £ ^ [28, 29]. By taking the scalar, and vector potentials as the 
Hulthen potential 



and using the Eq. (31), we get 



d 2 (f)(r) < 2 _ 2 2m (S - m,) + 2EV 2 mi So -m 2 + V 2 - S 2 
dr 2 + I m ° + e'A-o - 1 + ( e r/ro - 1)2 

+ l)e r / ro 



r 2( e r/r _ !)2 

By using a new variable e _r//r ° = s, Eq. (33) becomes 



r) = (31) 
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d 2 (j>(s) + 1-s d<t>{s) | j r 2 (E 2 -m 2 ) + 2r 2 (m (,g - mQ + EV ) 
ds 2 s(l — s) rfs | s 2 s(l — s) 

rgCm^So-mQ + ^-go 2 ) 1(1 + 1) ] 
+ (13^ s( i_ s)2 j^- U ' ^ 

By using the new parameters 



a(mi) = r)(mi)r , 

r) 2 (rrii) = (mi — m o) 2 — E 2 , 

&Vi) = r 2 [2EV -2S (m 1 -m )], 

(3 2 2 { mi ) = r 2 [2EV - 2m (m 1 - S )\ , 

v 2 {m x ) = -a 2 + a 2 (m 1 )+/3 2 (m 1 )-/3 2 (m 1 ) + z/ 2 , (33) 

where v(mi)(m\ — > 0) = v ,r](m 1 )(m 1 — > 0) = 77, and a = r)r , and comparing Eq. (34) with 
Eq. (11), we get the following parameter set given in Section II 



Oil = 1 , 
a 2 — 1 , 

"3 = 1, 
Ct4 = , 

"6 = 6 + 1) 
"8 = 6, 

aio = 1 + 2^ 

"12 = , 



£i = « 2 (mi) + ,3 2 (mi) + z/ 2 (mi) 
6 = 2a 2 + /? 2 (m 1 )-£(£+l) 
6 = « 2 

0^1 = -6 

a 9 = 6 - 6 + 6 + i 



(34) 



2 + 2(^1-6 + 6 + 1 + Via) 

«i3 = -1 - ( v^i-6 + 6 + i + v 7 ^) 

where z/ 2 = — V 2 ) , and r} 2 = m\ — E 2 in the above equations. 

We can easily get the energy eigenvalue equation of the Hulthen potential by using Eq. 
(20) 



/3 2 (m 1 )-£(£+l)-n 2 -(2n-l)^ 

a= 2(^TV) ' (35) 

where 5' — \ + \\j(2t + l) 2 + 4z/ 2 (mi) . We list some energy eigenvalues in Table I, and 
Table II for the case of constant mass, and the one of spatially dependent mass, respectively. 



To compare our results, we have used the values of the parameters given in Ref. [32]. E a 
denotes the energy eigenvalues of the particle, and E p denotes the one of the antiparticle in 
Table I, and Table II. 

According the result obtained in Eq. (37), we give easily the eigenvalue equation in the 
case of constant mass 



g(mi = 0) - £(£ + 1) - n 2 - (2n - 1)5 \ mi = 0) 
a= 2(n + *'(mi = 0)) ' (36) 

which is the same with the result obtained in Ref. [28]. 

The corresponding eigenfunctions of the Hulthen potential is written by using Eq. (26), 
and Eq. (35) 



4>[r) = A n e~ ar/ra (1 - e - r/ro ) 1+s ' P r {2a ' 1+2S '\l - 2e~ r/r °) , (37) 

where A n is a normalization constant. 

Finally, the eigenfunctions in the case of constant mass are written by using Eq. (39) 



0(r) = A' n e- ar/ro (1 - e - r/ro ) 1+s " Pi 2a ' 1+2S "\l - 2e- r/ro ) . (38) 

where 5" = \+\sJ(2e + l) 2 + ^v 2 {m x = 0) . The Jacobi polynomials P^ 2a ' 1+26 "\l-2e~ r ' r °) 
in the last result can be written in terms of hypergeometric function 2 i ? i(— n, n + 2a + 25 ' + 
2, 2a; s) , which gives the same result obtained in Ref. [28]. 

The normalization constant in Eq. (39) is obtained from the the normalization condition 

/ \<P(r)\ 2 dr = 1 , (39) 
Jo 

By introducing a new variable as x = 1 — 2e~ r / r °, we have from Eq. (41) 

Kl'^b j^^-xf a -\l + x){l + xfP^^\x)P^^\x)^l, (40) 
where j3 — 1 + 25' . By using the following identities [30, 31] 
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2n(C + (' + n)(( + C + 2n- 2)P«-^(x) = (C + C + 2n - 1)(( 2 - C*)j£i C) (*) 
+ (C + C' + 2n - 1)(C + C + 2n)(C + C' + 2n - 2)xP^P(x) 

-2(C + n - 1)(C + n - 1)(C + C + 2n)p£' 2 C,) (^) , (41) 

and 



wc+c 

we obtain the normalization constant 



n!a- 



(2a + /3 + 2n + 2)(2a + /3 + 2n) 



r(2a + /5 + n+ 1) 



.(43) 



n ^\ ' 4n(n + l + 2a + /3) + 2(l + /?)(2a + /3)r(2a + n+l)r(/5 + n+l) 

By following the same procedure, the normalization constant A' n in the eigenfunctions of 
the case of constant mass is obtained as A' n = A n (f3 — > 1 + 25") in Eq. (43). 



IV. CONCLUSION 



We have approximately solved the Klein-Gordon equation for the Hulthen potential for 
any angular momentum quantum number in the position-dependent mass background. We 
have found the eigenvalue equation, and corresponding wave functions in terms of Jacobi 
polynomials by using NU-method within the framework of an approximation to the centrifu- 
gal potential term. We have also obtained the energy eigenvalue equation, and corresponding 
eigenfunctions for the case of the constant mass. Results for the case of constant mass are 
the same with the ones obtained in Ref. [28]. 
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TABLE I: The energy eigenvalues of vector, and scalar Hulthen potential for m = 1 , and m\ = 0. 



V = S = 1 





n 
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E a 


771 (i 

tip 


T7> b 


7—1 f) 


7— ' r 

E a C 


ftp 




1 







-0.6000000 


1.0000000 


-0.6000000 


1.0000000 


-0.6000000 


1.0000000 




1 

= S 


= 2 


1 
















1 







r\ '~rr\'~r'\ n/~»o 

-0.7071068 


0.7071068 


-0.7071068 


0.7071068 


-0. 7071068 


n ^7r\ 77 1 r\/^o 

0.7071068 




-i 
l 




1 


n 01 A cia n7 


n aqji 7-1 /| 

U.9o41 1 14 












2 
2 




u 

1 


-U.214y4U/ 


U.9o41 / 14 


n 11 /in/i i n 

-U.zl4y4iu 


n no A 1 71 n 
0.9841 / 10 


n 01 /in/1 1 n 
-U.214y41U 


n no /i 1 71 n 
0.9841 / 10 


V 


2 

= S 


= 3 


2 


— 


— 


— 


— 


— 


— 




1 







-0.7637079 


0.3021695 


-0.7637080 


0.3021690 


-0.7637080 


0.3021690 




1 




1 


-0.4114378 


0.9114378 












2 




n 

u 


A /| 1 1 /| OTO 

-0.4114378 


0.9114378 


A /I 1 1 A oon 

-{J.4114ao(J 


n m i a oon 

0.9114380 


r\ /ill /loon 
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n ni 1 i oon 

0.9114380 
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0.6000000 


0.6000000 
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0.6000000 
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n n c\c\c\c\c\c\ 

0.6000000 
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0.2844158 


0.9942727 


0.284416 


0.994273 


0.2844160 


0.9942730 



°our results 

^results obtained in Ref. [32] 

"results obtained in Ref. [33], and Ref. [34] 
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TABLE II: The energy eigenvalues of vector, and scalar Hulthen potential for rri\ / 0. 
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continued 
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